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Using a data sample of 10.5 fb 1 collected by the Belle detector, three-body charmless decays _B+ — > K + h + h have been studied. 
The following branching fractions have been obtained: B{B+ -» K"+7T~7T+) = (64.8 ± 10.0± 7.0) X 10 -6 and B(B+ -+ = 
(36.5 ± 6.1 ± 5.5) x 10" 6 . The upper limits for other combinations of charged kaons and pions have been placed. Analysis of the 
intermediate two-body states gives evidence for production of scalar resonances in charmless B decays. 



1 Introduction 

Until recently charmless B decays were mainly studied 
via their two-body decay modes because of the large com- 
binatorial background in multibody channels. In this 
work we attempt to study B + — > K + h + h~ decays (h 
stands for a charged pion or kaon) without any assump- 
tion about the intermediate hadronic resonances. Refer- 
ence to the charge conjugate states is implicit throughout 
this paper unless explicitly stated otherwise. 

Analysis i^ based on the data set collected by the 
Belle detector El at KEKB □, the asymmetric B-factory 
at KEK. It consists of 10.5 fb" 1 taken at the T(4S) 
resonance and 0.6 fb -1 taken below the BB production 
threshold for continuum studies. 

2 Event selection 

Charged tracks are required to satisfy a set of track qual- 
ity cuts based on the the average hit residual and the 
impact parameters in both the r — <fi and r — z planes. 
We require that the transverse track momenta be greater 
than 100 MeV/c to reduce low momentum combinatorial 
background. 

The candidate events are identified by using the 
beam-constrained mass Mt c = sj (y/s/2) 2 — P B 2 and the 
calculated energy difference AE = E* B — s/s/2 where E* B 
and P B are the energy and 3-momentum of the B candi- 
date in the T(4S) rest frame. We define the B signal re- 
gion as: 5.271 < M bc < 5.289 GeV/c 2 ; \AE\ < 40 MeV. 

The most important issues of this analysis are ef- 
ficient kaon identification over the whole momentum re- 
gion and the suppression of the large combinatorial back- 
ground which is dominated by qq continuum production. 

Since two B mesons produced from the T(4«S) decay 
are nearly at rest in the CMS frame, the angles of the 
decay products of two B's are uncorrelated and the event 
looks spherical. In contrast, hadrons from continuum qq 



events tend to exhibit a two-jet structure. 

We calculate the angle, Qtkt , between the thrust axis 
of the B candidate and that of the rest of the event. The 
distribution of | cos(0Thr)\ is strongly peaked near 1.0 for 
qq events while it is nearly flat for BB events. We require 
I cos(9xhr)\ < 0.80 for all modes under consideration. 

We also form a Fisher discriminant u with the mo- 
mentum scalar sum of charged particles and photons in 
nine cones of increasing polar angle around the thrust 
axis of the B candidate and the angle of the thrust axis 
of the candidate with respect to the beam axis. We com- 
bine the Fisher discriminant and cos (9b), where Ob is 
the angle between the B candidate momentum and the 
beam axis, in a single variable by taking a product of the 
corresponding probability density functions and impose 
a cut on the likelihood ratio. 

Separation of kaons and pions is accomplished by 
combining the responses of the ACC and the TOF with 
dE/dx measurements in the CDC. The combined re- 
sponse of the three systems provides high quality K/tt 
separation in the laboratory momentum range up to 3.5 
GeV/c. At large momenta only the ACC and dE/dx are 
used since here TOF provides no significant separation 
of kaons and pions. 

3 Analysis 

In this analysis at least one charged track was required to 
be positively identified as a kaon. Then we considered all 
possible combinations of charged kaons and pions. The 
analysis of K + tt + tt" and K + K + K~ final states is de- 
scribed in detail in two following subsections. The re- 
sults on the study of other combinations are summarized 
in Table § 



3.1 £r 



To study this final state we require one track to be 
positively identified as a kaon and two tracks consis- 
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Figure 1: The M bc (top) and AE (bottom) distributions for se- 
lected -B+ — > K + tt + tt~ candidates. Points are data and his- 
tograms are Monte Carlo expectation for background. Curves are 
fit to the data. 



Figure 2: The K + ir~ (top) and 7T+7T - (bottom) invariant mass 
spectra for candidates from the B signal box. Histograms are data 
and curves are fitting results. The hatched regions are the back- 
ground estimates from sidebands. 



tent with the pion hypothesis. The contributions from 
B+ -► D°ir+ where D° -► K + tt~ and B+ -► J/ip(ip')K+ 
where J/ip(ip') — ► submodes were excluded from 

the analysis by imposing cuts on the invariant masses 
of two intermediate particles: \M(K + ir~) — 1 .865 1 > 
0.100 GeV/c 2 ; \M(h + h~) - 3.097| > 0.070 GeV/c 2 ; 
\M(h + h~) - 3.686| > 0.050 GeV/c 2 where h+ and h~ 
are pion candidates. The latter two submodes contribute 
due to the muon-pion misidentification and in these cases 
we use the muon mass for the M(h + h~) calculation. 

The Mb c and AE distributions after excluding these 
signals are presented in Fig. |l|. A significant enhancement 
in the B signal region can be seen in both distributions. 

To determine the intermediate states which con- 
tribute to the observed signal, we plot the K + n~ and 
7r + 7r~ invariant mass spectra as shown in Fig. |^. The 
dashed regions in Fig. || show the background spectra 
determined from Mt, c and AE sidebands. We fit the sig- 
nal in the K + n~ invariant mass spectrum to the non- 
coherent sum of two relativistic Breit-Wigner functions. 
The parameters of the first Breit-Wigner were fixed to be 
equal to those of the if* (892) meson and the parame- 
ters of the second one referred to as if^-°(1400) were free 
during the fit. The signal in the 7r + 7r~ invariant mass 
spectrum is fitted to the non-coherent sum of three rel- 
ativistic Breit-Wigner functions. The parameters of the 
first and second Breit-Wigners were fixed to be equal to 
those of p°(770) and / O (980) mesons respectively and the 
parameters of the third one referred to as /x(1300) were 
free during the fit. Results of the fit are summarized in 
Table g. 



3.2 B+ -> K+K+K- 

We select B candidates formed from three charged tracks 
identified as kaons. The contribution from the Cabibbo 
suppressed B + — ► D°K + decay where D° — > K + K~, 
is excluded from the analysis by imposing cuts on the 
K+K- invariant mass: \M{K + K~) - 1.865| > 0.025 
GeV/c 2 . To suppress the background from tt— K misiden- 
tification we reject all candidates if the invariant mass of 
any two oppositely charged tracks is consistent within 
2a with the D — > Ktt hypothesis with the kaon hypoth- 
esis assigned to be either of the two tracks. No particle 
identification information is used in this veto. 

The resulting Mt, c and AE distributions are pre- 
sented in Fig. ^|. A large enhancement in the B signal 
region can be seen in both distributions. 

To determine the intermediate states which con- 
tribute to the observed signal, we examine the K + K~ 
invariant mass spectrum which is shown in Fig. [|. The 
dashed region in Fig. |] shows the background spectrum 
determined from Mb c and AE sidebands. 

The K + K~ invariant mass spectrum is character- 
ized by two distinct structures. The first one around 
1.0 GeV/c 2 is very narrow and corresponds to the 
0(1020) meson. The other one around 1.5 GeV/c 2 which 
we refer to as /x(1500) is very broad. We fit the signal 
in the K + K~ invariant mass spectrum to the relativis- 
tic Breit-Wigner function. The K + K~ invariant mass 
region correspondig to the 0(1020) meson was excluded 
from the fit. The results of the fit are summarized in 
Table @. 
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Figure 3: The M bc (top) and AE (bottom) distributions for se- 
lected B + — > K + K + K~ candidates. Points are data and his- 
tograms are Monte Carlo expectation for background. Curves are 
fit to the data. 

4 Results 

The final signal yield for each submode was extracted 
from the fit to the beam constrained mass distributions 
after cuts on the invariant mass of two intermediate par- 
ticles are applied. We fit the M(, c distributions to the 
sumjpf a Gaussian function for signal and ARGUS func- 
tion □ for background. We assume that yield in each bin 
of K + ir~ /it + it~ is entirely due to a single resonance. 

To reduce the systematic error in the branching 
fractions we normalize our results to the B + — ► D°ir, 
D — ► K + tt~ signal. Similarly, because of poor knowl- 
edge of the absolute branching fractions and the uncer- 
tainty in the interpretation of some of intermediate reso- 
nances, we calculate only the products of the branching 
fractions rather than their absolute values: 

N x 1 

&B+~,Xh+ X Bx^h+h- = T7- X Bb^Dtt X Bd^Ktt X - 

N Dv d 

where X denotes a h + h~ resonance state, N x and 
are the number of observed events for a mode under 
study and for the reference process respectively, Bb^dtt 
and Bd^k-it are branching fractions for B + —> D°ir and 
D° — > K + n~ respectivelycl and S is the overall correction 
factor. 

5 Discussion & Conclusions 

The high quality of tt/K separation allowed us to mea- 
sure for the first time the branching ratios of the three- 
body decays B+ ->■ K+tt+tt' and B+ -> K + K + K~ 
without any assumptions about particular intermedi- 
ate mechanisms. CLEOQ and BaBar 13 have previously 



Figure 4: The K + K~ invariant mass spectrum for candidates from 
the B signal box. The hatched region is the background estimates 
from sidebands. 



Table 1: Results of the fit to the K+tt , 7r + 7r and K + K in- 
variant mass spectra. 



State 


M, GeV/c 2 


T, GeV 


Yield 


if*°(892) 


0.896 


0.051 


8.7 ±4.7 


#|?(1400) 


1.31 ±0.05 


23+° 10 


28.3 ±7.4 


p°(770) 


0.769 


0.151 


10.2 ±5.8 


/ O (980) 


0.980 


0.060 


22.8 ±4.9 


/x(1300) 


1.33 ±0.03 


13+ 007 

u - io -0.04 


17.7 ±6.6 


0(1020) 


1019.4 




14.1 ±4.6 


/x(1500) 


1.43 ±0.08 


41+° 18 


32.8 ±7.6 



placed upper limits on non- resonance three-body decays. 
The reported numbers for B + — > K + -k + ii~ (CLEO: 
< 28 x 10~ 6 , BaBar: < 66 x 10~ 6 ) and B+ -> K+K+R- 
(CLEO: < 38 x 10~ 6 ) are considerably lower than those 
presented in this paper. Comparison of the applied se- 
lection criteria shows that CLEO and BaBar restricted 
their analysis to the region of the invariant masses above 
2 GeV/c 2 for any pair of the particles. Assuming the 
phase space distribution of the invariant masses, they 
obtained the limits quoted above. Similar analysis in our 
case gives consistent results. 

The upper limits obtained for the if~7r + 7r + , 
K + K~n + and K + K + tt~ modes, are considerably better 
than previous results by CLEOH and OPALEl. A search 
for the K + K + ir~ mode is of particular interest since in 
some extensions of the Standard Model this branching 
fraction is predicted to be significantly enhanced. 

A clear signal has been found in the channel B + — ► 
i^°(1400)7r+, 2f£>(1400) -> K+tt-. Among the reso- 
nances in this region, ifi(1270), if i (1400), if*(1410), 
if* (1430) and if|(143Q), the first two do not decay into 
Kir at all. The branching fraction of if* (1410) — > Kit is 
too small to make its contribution essential. Thus, only 
two candidates remain: if* (1430) and if* (1430). How- 
ever, in the factorization approximation which seems ad- 
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Table 2: Summary of the results on the search for B + —* K+h+h 
decays. 



Table 3: Summary of the results on the search for intermediate 
resonances in B+ — > if +7r + 7r — and B + — > K + K~ states. 













Mode 


S 


Yield 


Sig., a 


Bxh x Z?/^ 


Mode 


5 


Yield 


Sig., a 


s(io- 6 ) 






0.94 


83.0±12.5 


8.6 


64.8±10.±7.0 




K*°(892)n+ 


0.92 


7-5+fi 


2.8 


< 11.6 


K+K+K- 


0.72 


48.0 ±8.0 


9.3 


36.5±6.1±5.5 




^°(1400)7r+ 


0.89 


29.11S5 


5.2 


23.8lg;i±5.8 


K~1T + TT + 


0.94 


7.8±|;| 




< 12.8 




p°(770)K+ 


0.67 


5 5+ 4 - 1 

°-°-3.4 


1.6 


< 13.5 


K+K-TT+ 


0.83 


9.71^ 




< 17.0 




f o (9S0)K+ 


0.79 




4.9 


16.4±^±2.8 


K+K+n- 


0.83 






< 5.2 




f x (1300)K+ 


0.79 


20 4+ 5 - 7 


4.8 


18.8±|-J±4.3 












fx{l500)K+ 


0.45 




5.5 


25.1±^±4.3 



equate in this case, the local production of the tensor 
if|(1430) is strictly zeroB. Therefore, its production is 
only possible due to corrections to the factorizable con- 
tribution and most probably is highly suppressed. 

In contrast, the scalar meson if* (1430) can be easily 
produced via the main factoxizable contribution and ac- 
cording to a recent estimated the branching fraction for 
B+ if*(1430)7r+ is even larger than that for B+ -> 
K°tt+. The prediction gives B(B+ -> K * (1430)tt+) x 
B(K*(U30) K+tt-) ~ 28 x 10~ 6 which is in good 
agreement with our result (23.8±l* ± 5 - 8 ) x 10 ~ 6 - 

A significant signal was observed for the first time 
in the decay mode B + — ► f o (980)K + with a prod- 
uct of branching fraction of B(B+ -> f o (980)K+) x 
S(/ O (980) -> tt+ti-) = (16.4±J;| ± 2.8) x 10~ 6 also indi- 
cates strong coupling of penguins with scalars. 

It is not easy to identify a broad state with a mass 
around 1500 MeV observed in the decay mode B + — ► 
K + K + K~ . If an attempt is made to describe the ob- 
served spectrum by a single resonance in addition to the 
0(1020), then it is natural to assume that it is a scalar 
ss state. Such a state is hardly compatible with either 
one of the scalars /Jjl370) and / o (1500) suggested by 
Particle Data Group a and having small probabilities of 
the decay into K K final state. Existence of a state with 
required properties is predicted in some potential mod- 
els Ed and evidenceJbr such a state was reported by the 
LASS experiment El. 

Similar uncertainties arise in the interpretation of the 
peak with a tt + tt~ mass about 1300 MeV in the Kirir 
system. Two candidates for such a state - /2(1270) and 
/ O (1370) exist. The assumption of the / O (1370) mech- 
anism with its rather small coupling El to 7r + 7r~ would 
lead to an unnaturally large branching ratio of the B de- 
cay. The parameters of the bump obtained from the fit 
are close to those of the /2(1270). However, as recently 
shown til, factorization leads to a very small branching 
ratio for B + — > /2(1270)iT + . If our observation is shown 



to be caused by the /2(1270), then it provides evidence 
for the nonfactorized contributions which could also ac- 
count for the large yield of rj'K in B decays. 

Any unambiguous conclusion about the substructure 
of the observed signal will require a larger data sample 
which will then allow us to perform the angular analysis 
of the decay products of the intermediate resonances. 
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